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Abstract Mineral malnutrition stemming from undiversified plant-based diets is a top global 
challenge. In C3 plants (e.g., rice, wheat), elevated concentrations of atmospheric carbon dioxide 
(eC02) reduce protein and nitrogen concentrations, and can increase the total non-structural 
carbohydrates (TNC; mainly starch, sugars). However, contradictory findings have obscured the effect 
of eC02 on the ionome — the mineral and trace-element composition — of plants. Consequently, 
C02-induced shifts in plant quality have been ignored in the estimation of the impact of global change 
on humans. This study shows that eC02 reduces the overall mineral concentrations (-8%, 95% 
confidence interval: -9.1 to -6.9, p<0. 00001) and increases TNC:minerals > carbon:minerals in C3 
plants. The meta-analysis of 7761 observations, including 2264 observations at state of the art FACE 
centers, covers 130 species/cultivars. The attained statistical power reveals that the shift is systemic 
and global. Its potential to exacerbate the prevalence of 'hidden hunger' and obesity is discussed. 
DOI: 10.7554/eLife.02245.001 
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Introduction 

Mankind's ultimate source of carbohydrates is atmospheric carbon dioxide (CO2) converted by photo- 
synthesis to sugars. The bulk of the terrestrial conversion of C02-to-carbohydrates is done by C3 
plants, which account for over three quarters of global primary production and for over 90% of Earth's 
plant species {Still and Berry, 2003). {If not stated otherwise, hereafter, terms 'plant(s)' and 'crop(s)' 
refer to C3 species). When exposed to CO2 concentrations twice the preindustrial level of ~280 ppm, 
plants increase the synthesis of carbohydrates by 19-46% {Leakey et al., 2009). Currently, CO2 
concentrations are reaching 400 ppm — the highest level since the dawn of agriculture and likely to be 
the highest since the rise of modern humans {Siegenthaler et al., 200S}. Within a single human lifespan, 
CO2 levels are projected to reach 421-936 ppm {IPCC, 2013). Will rising CO2 concentrations — one 
of the most certain and pervasive aspects of global climate change — alter the quality of crops and wild 
plants? Will the C02-induced stimulation of carbohydrate synthesis increase the carbohydrates-to- 
minerals ratio in crops? Can such shifts in crop quality affect human nutrition and health? 

Elevated CO2 effects on plant quantity (productivity and total biomass) have been extensively 
studied and show higher agricultural yields for crops, including wheat, rice, barley, and potato. But 
eC02 effects on plant quality, and possible cascading effects on human nutrition, have been largely 
ignored in the estimation of the impact of eC02 on humans. Notably, IPCC (2007, 2013) and AAAC 
Climate Science Panel (2014) include direct CO2 effects (e.g., ocean acidification) in their climate 
change assessments but do not mention any CO2 effects on crop or wild plant quality. However, it is 
unwarranted to assume that plants will balance the increased carbohydrate synthesis with other adjust- 
ments to their physiology to maintain the nutritional quality for their consumers in a state of unper- 
turbed homeostasis. The stoichiometry — the relative ratios of chemical elements — in plants is plastic 
and, to a considerable degree, reflects their environment {Sterner and Elser, 2002). However, detecting 
C02-induced shifts in plant quality is challenging for several reasons. First, plant quality involves 
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eLlfe digest Rice and wheat provide two out every five calories that humans consume. 
Like other plants, crop plants convert carbon dioxide (or CO2) from the air into sugars and 
other carbohydrates. They also take up minerals and other nutrients from the soil. 

The increase in CO2 in the atmosphere that has happened since the Industrial Revolution is 
thought to have increased the production of sugars and other carbohydrates in plants by up to 
46%. CO2 levels are expected to rise even further in the coming decades; and higher levels of 
CO2 are known to lead to lower levels of proteins in plants. But less is known about the effects 
of CO2 levels on the concentrations of minerals and other nutrients in plants. 

Loladze has investigated the effect of rising CO2 levels on the nutrient levels in food plants by 
analyzing data on 130 varieties of plants: his dataset includes the results of 7761 observations made 
over the last 30 years, by researchers around the world. Elevated CO2 levels were found to reduce 
the overall concentration of 25 important minerals — including calcium, potassium, zinc, and iron — in 
plants by 8% on average. Furthermore, Loladze found that an increased exposure to CO2 also 
increased the ratio of carbohydrates to minerals in these plants. 

This reduction in the nutritional value of plants could have profound impacts on human health: 
a diet that is deficient in minerals and other nutrients can cause malnutrition, even if a person 
consumes enough calories. This type of malnutrition is common around the world because many 
people eat only a limited number of staple crops, and do not eat enough foods that are rich in 
minerals, such as fruits, vegetables, dairy and meats. Diets that are poor in minerals (in particular, 
zinc and iron) lead to reduced growth in childhood, to a reduced ability to fight off infections, and 
to higher rates of maternal and child deaths. 

Loladze argues that these changes might contribute to the rise in obesity, as people eat 
increasingly starchy plant-based foods, and eat more to compensate for the lower mineral levels 
found in crops. Looking to the future, these findings highlight the importance of breeding food 
crops to be more nutritious as the world's CO2 levels continue to rise. 
DOI: 10.7554/eLife.02245.002 



multiple nutritional currencies, for example macronutrients (carbohydrates, protein, and fat) and 
micronutrients (minerals, vitamins and phytonutrients). Assessing relative changes within and 
among multiple currencies requires significantly more effort and funding than measuring only 
plant quantity (e.g., yield). Second, plant quality, including the plant ionome — all the minerals and 
trace-elements found in a plant {Lahner et a/., 2003; Salt et al., 2008) — is inherently variable; and 
measurement imprecisions further amplify the variability. For example, Stefan et al. (1997) report 
the accuracy test for 39 facilities that analyzed samples of the same plant tissues: the inter-laboratory 
variance was 6.5% for N, but twice as large for phosphorus (P) and calcium (Ca), and reached 1 30% 
for sodium (Na). Therefore, C02-induced changes in the plant ionome (the signal) can be easily 
lost amid highly variable data (the noise), especially when such data are limited and sample sizes 
are small. However, it is important to bear in mind that a low signal-to-noise ratio does not imply 
that the signal is practically insignificant, especially if it is global and sustained — a point revisited 
in the 'Discussion'. 

Elusive CO2 effect on the plant ionome: contradictory findings 

The first empirical evidence of lower mineral content in plants exposed to eC02 appeared at least over 
a quarter century ago (e.g.. Porter and Grodzinski, 1984; Peet et al., 1986; O'Neill et al., 1987). 

Physiological mechanisms responsible for the overall decline of plant mineral content — with expected 
changes being non-uniform across minerals — have been proposed: the increased carbohydrate 
production combined with other eC02 effects such as reduced transpiration {Loladze, 2002; McGrath 
and Lobell, 2013). However, most of the experimental evidence showing C02-induced mineral 
declines came from artificial facilities, mainly closed chambers and glasshouses, and many results were 
statistically non-significant. This led some research groups to challenge altogether the notion of lower 
mineral content in plants exposed to eC02 in field conditions. Such conditions are most accurately 
represented in Free-Air Carbon dioxide Enrichment (FACE) centers, which have been established in at 
least 1 1 countries. 
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In the grains of rice harvested at four FACE paddies in Japan, Lieffering et al. (2004) found no 
decline in any of the minerals but lower N content. The result disagreed with Seneweera and Conroy 
(1997), who were the first to report lower iron (Fe) and zinc (Zn) in grains of rice grown at eCOa and 
warned that altered rice quality can negatively affect developing countries. Lieffering et al. (2004), 
however, argued that the result of Seneweera and Conroy (1997) could be an artifact of growing 
rice in pots, which restrict rooting volumes. They hypothesized that in FACE studies, which provide 
unrestricted rooting volumes, plants would increase uptake of all minerals to balance the increased 
carbohydrate production. This hypothesis, however, found no support in the FACE studies of Pang et al. 
(2005) and Yang et al. (2007) (carried out in China and latitudinally not very far from the study in 
Japan), who found that eC02 significantly altered the content of several minerals in rice grains. 

The contradictory results coming from these studies on rice seem perplexing, especially in light of 
the very robust effect that eCOa has on N in non-leguminous plants. Elevated CO2 reduces N concen- 
trations by 10-18% systemically throughout various tissues: leaves, stems, roots, tubers, reproductive 
and edible parts, including seeds and grains {Cotrufo et al., 1998; Jablonski et al., 2002; Taub et al., 
2008). If the increased carbohydrate production dilutes the nutrient content in plants, why does the 
dichotomy seem to exist between the responses of N and minerals to eC02? In addition to the carbo- 
hydrate dilution and reduced transpiration, eC02 can further lower N concentrations in plants by: (1) 
reducing concentrations of Rubisco — one of the most abundant proteins on Earth that comprises a 
sizable N-pool in plants {Drake et al., 1997), and (2) inhibiting nitrate assimilation {Bloom et al., 
2010). Hence, it is reasonable to expect the effect of eC02 on N to be larger and, thus, easier to dis- 
cern than its effect on most minerals. The stronger signal for N, combined with the plentiful and less 
noisy data on this element, can help explain why by the end of last century the effect of eC02 on N 
had been already elucidated {Cotrufo et al., 1998), but its effect on minerals has remained elusive. 

The obscure nature of the effect of eC02 on minerals becomes particularly apparent in the largest 
to date meta-analysis on the issue by Duval et al. (2011), who fragmented data from 56 eC02 studies 
into 67 cases. In 47 of the cases, the effect of eC02 on minerals was statistically non-significant, that is 
the 95% Confidence Interval (CI) for the effect size overlapped with 0. The remaining 20 cases were 
statistically significant but showed no pattern: for example, Fe increased in grasses but decreased in 
trees, Zn increased in roots but decreased in stems, while in grains only sulfur (S) decreased. Duval 
et al. (2011) concluded: "A major finding of this synthesis is the lack of effect of CO2 on crop grains 
nutrient concentration". This would imply laying to rest the hypothesis that eC02 consistently alters 
the plant ionome and would render mitigation efforts to combat declining crop mineral concentra- 
tions in the rising CO2 world unnecessary. However, a closer examination of the results of Duval et al. 
(2011) reveals that every statistically significant increase in mineral concentrations was obtained by 
bootstrapping a sample of size 2, 4 or 5 — a recipe for generating invalid 95% CIs. loannidis (2005) 
showed that false research findings, stemming from small sample sizes and associated low statistical 
power, are a persistent problem in biomedical sciences. 

'Power failure' and the plant ionome 

Calling the problem as 'power failure', Button et al. (2013) emphasized that the probability of a 
research finding to reflect a true effect drops drastically if the statistical power is reduced from 0.80 
(considered as appropriate) to low levels, for example <0.30. Since the power of a statistical test drops 
non-linearly with the effect size, a sample size that is sufficient for detecting a 15% effect, for example 
a decline in N content, can be inadequate for detecting a 5% effect, for example a decline in a mineral 
content. Considering that the standard deviation of mineral concentrations in a plant tissue can reach 
25% {Duquesnay et al., 2000; Lahner et al., 2003), the 5% effect size standardized as Cohen's d is 
d = 5/25 = 0.2. A t test applied for d = 0.2 to a sample size of 3-5 — a typical size used in eC02 
studies — yields the power of 0.06-0.1 0 {Paul et al., 2007). (Unfortunately, MetaWin {Rosenberg et al., 
2000), a statistical package routinely used in meta-analytic and other CO2 studies in ecology, provides 
neither a priori nor post-hoc power estimates.) Such a small power not only raises the probability of 
obtaining a false negative to 90-94% but also increases the likelihood that a statistically significant 
result does not reflect a true effect {Button et al., 2013). 

Answering questions with adequate power 

As of this writing, researchers on four continents have generated data sufficient for answering with an 
adequate statistical power, the following questions: 
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1. Does eC02 shift the plant ionome? If yes, what are the direction and magnitude of shifts for indi- 
vidual chemical elements? How does the effect of eC02 on N compares to its effect on minerals? 

2. Do FACE studies differ principally from non-FACE studies in their effect on the plant ionome? 

3. Do the plant ionomes in temperate and subtropical/tropical regions differ in their response to eCOa? 

4. Do the ionomes of photosynthetic tissues and edible parts differ in their response to eCOa? How 
does eC02 affect the ionomes of various plant groups (woody/herbaceous, wild/crops, C3/C4) and 
grains of the world's top C3 cereals — wheat, rice, and barley? 

Results 

For brevity, hereafter 'minerals' refer to all elements except C, hydrogen (H), oxygen (O), and N. 
All results are for C3 plants except when noted otherwise. 

Power analysis 

Plotting the effect sizes (with 95% CIs) for the 25 minerals against their respective statistical power 
reveals a clear pattern {Figure 1). In the very low power (<0.20) region, the noise completely hides 
the COa-induced shift of the plant ionome. In the low power region (<0.40), the shift still remains 
obscure. As the statistical power increases, so does the likelihood that a statistically significant 
result reflects true effect and, consequently, the direction and the magnitude of the CO2 effect on 
minerals become increasingly visible in the higher power regions of the plot. 

To increase the likelihood of reporting true effects, only results with the statistical power >0.40 are 
reported in this section. However, Figure 1-source data 1 lists all the results together with their 
p-values irrespective of the statistical power (e.g., results for chromium (Cr) or the bean ionome are not 
shown here due to low power, but are listed in Figure 1-source data 1). 




statistical significance: 

• Significant increase 

• Nonsignificant ctiange 

• Signi^cant decrease 



-10 0 10 

Effect size (%) 



Figure 1 . Statistical power and the effect of CO2 on the plant ionome. The effect of elevated atmospheric CO2 concentrations (eC02) on the mean 
concentration of minerals in plants plotted (with the respective 95% confidence intervals [CI]) against the power of statistical analysis. The figure reflects 
data on 25 minerals in edible and foliar tissues of 125 C3 plant species and cultivars. The true CO2 effect is hidden in the very low and the low power 
regions. As the statistical power increases, the true effect becomes progressively clearer: the systemic shift of the plant ionome. 
DO!: 10.7554/eLife.02245.003 

The following source data are available for figure 1 : 
Source data 1. Supportive data for Figures 1-8. 

DO!: 10.7554/eLife.02245.015 



Loladze. eLife 2014;3:e02245. DOI: 10.7554/eLife.02245 



4 of 29 



Research article 



Ecology | Epidemiology and global health 



CO2 effect on individual elements 

Across all the data, eC02 reduced concentrations of P, potassium (K), Ca, S, magnesium (Mg), Fe, Zn, 
and copper (Cu) by 6.5-10% (p<0.0001) as shown on Figure 2. Across all the 25 minerals, the mean 
change was (-8%, -9.1 to -6.9, p<0. 00001). Only manganese (Mn) showed no significant change. 
It is not clear whether the oxygen-evolving complex (OEC) demands for Mn separate this mineral 
from the pattern of declines exhibited by other minerals. Among all the measured elements, only C 
increased (6%, 2.6 to 10.4, p<0.01 ). The sharp difference between the responses of C and minerals 
to eC02 is expected if a higher carbohydrate content drives the change in the plant ionome: for most 
plant tissues, the dilution by carbohydrates lowers the content of minerals while having little effect 
on C [Loladze, 2002). (This also suggests that the increase in C concentrations found here could 
be a result of a higher content of lipids or lignin — the two sizable plant compounds that are very 
C-rich [-60-75% C].) 

The patterns of change within edible and foliar tissues are similar: N, P, Ca, Mg, Zn, and Cu declined 
significantly in both tissues {Figures 3, 4). Aside from Mn, only K showed no significant decline in the 
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Figure 2. The effect of CO2 on individual chemical elements in plants. Change (%) in the mean concentration of 
chemical elements in plants grown in eC02 relative to those grown at ambient levels. Unless noted otherwise, all 
results in this and subsequent figures are for C3 plants. Average ambient and elevated CO2 levels across all the 
studies are 368 ppm and 689 ppm respectively. The results reflect the plant data (foliar and edible tissues, FACE 
and non-FACE studies) from four continents. Error bars represent the standard error of the mean (calculated using 
the number of mean observations for each element). The number of mean and total (with all the replicates) 
observations for each element is as follows: C(35/169), N(1 40/696), P(1 52/836), K(1 28/605), Ca(1 39/739), S(67/373), 
Mg(123/650), Fe(125/639), Zn(123/702), Cu(124/612), and Mn(1 01/493). An element is shown individually if the 
statistical power for a 5% effect size for the element is >0.40. The 'ionome' bar reflects all the data on 25 minerals 
(all the elements in the dataset except of C and N). All the data are available at Dryad depository and at GitHub. 
Copies of all the original sources for the data are available upon request. 
DOI: 10.7554/eLife.02245.004 
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Foliar tissues 
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Figure 3. The effect of CO2 on foliar tissues. Change (%) in the mean concentration of chemical elements in foliar 
tissues grown in eC02 relative to those grown at ambient levels. Average ambient and eC02 levels across all the 
foliar studies are 364 ppm and 699 ppm respectively. Error bars represent 95% CI. For each element, the number of 
independent mean observations, m, is shown with the respective statistical power For each plant group, m equals 
the sum of mean observations over all the minerals (not including C and N) for that group. Elements and plant 
groups for which the statistical power is >0.40 (for a 5% effect size) are shown. 
DOI: 10.7554/eLife.02245.005 



edible tissues (on Figure 1, it is visible as one of the only two black 95% CI in the 'High Power' region). 
In the foliar tissues, Mg declined the most (-12.3%, -16 to -8.7), which is congruent with the 
hypothesis of McGrath and Lobell (2013) that Mg should exhibit a larger decline in photosyn- 
thetic tissues because 'chlorophyll requires a large fraction of total plant Mg, and chlorophyll 
concentration is reduced by growth in elevated CO2'. However, the 95% CIs for Mg and for most 
other minerals overlap. A richer dataset would shed more light on the issue of Mg in photosyn- 
thetic tissues. 

As expected, among all elements N declined the most (-15%, -17.8 to -13.1, p<0.00001) {Figure 2), 
matching very closely previous findings [Figures 3-6): the 17-19% decline in leaves found by Cotrufo 
et al. (1998) and the 14% decline in seeds found by Jablonski et al. (2002). Since the contents of N 
and protein correlate strongly in plant tissues, the lower N in edible tissues {Figure 4) corroborates the 
protein declines in crops found by Taub et al. (2008). 

FACE vs non-FACE studies 

With respect to the types of experiments, the CO2 effect on the plant ionome is surprisingly robust: 
in both the FACE and the non-FACE studies eC02 significantly reduced N, P, K, Ca, S, Mg, and Zn 
{Figures 5, 6). The high cost of CO2 required for running free-air experiments led to a much lower 
average level of eC02 in the FACE studies (560 ppm) cf. 732 ppm in the non-FACE studies. It is plau- 
sible that the lower levels of CO2 in the FACE studies contributed to a smaller overall mineral decline 
(-6.1%, -7.8 to -4.4) cf. (-8.7%, -10.1 to -7.4) for the non-FACE studies. In both the FACE and the 
non-FACE studies, the overall mineral concentrations declined significantly in herbaceous plants and 
crops, foliar and edible tissues, including wheat and rice {Figures 5, 6). 

Geographical analysis 

The CO2 effect on the plant ionome appears to be pervasive throughout latitudes {Figures 7, 8). With the 
exception of three small centers (in Bangladesh, Japan, and the UK), the mean mineral concentrations 
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Figure 4. The effect of CO2 on edible tissues. Change (%) in the mean concentration of chemical elements in edible 
parts of crops grown in eC02 relative to those grown at ambient levels. Average ambient and elevated CO2 levels across 
all the crop edible studies are 373 ppm and 674 ppm respectively. Other details are in the legends for Figures 2 and 3. 
DOI: 10.7554/eLife.02245.006 



declined in every FACE and open top chamber (OTC) center on four continents. The mineral 
decline in the tropics and subtropics (-7.2%, -10.4 to -4.0, p<0.0001) is comparable to the decline in 
the temperate region (-6.4%, -7.9 to -5.0, p<0. 00001). A finer regional fragmentation currently is not 
possible due to lack of data for Africa, South America, Russia, and Canada. For many existing centers 
the data are limited and yield a low statistical power. 

Germany leads the world in the FACE and OTC data generation with the largest number of mean 
observations of mineral concentrations (285), followed by the USA (21 8) {Figure 8). Though Australia 
generated only 30 mean observations, it stands out in the exceptional precision of some of its studies: 
the wheat experiments of Fernando et a/. (2014) employed an unprecedented for FACE studies 48 
replicates (for this reason, the study is easily identifiable on Figure 9). 

CO2 effect on various plant groups and tissues 

Since eC02 does not stimulate carbohydrate production in C4 plants to a degree that it does in C3 plants, 
one would expect a milder CO2 effect on minerals for C4 plants. Indeed, no statistically significant effect 
was found on the ionome of C4 plants {Figure 8). Note, however, that the very limited data on this 
plant group are insufficient for deducing the absence of the effect; rather, it is likely that the effect 
size <5% for C4 plants. 

The CO2 effect on the C3 plant ionome shows its systemic character through the analysis of various 
plant groups and tissues {Figures 3, 4 and 8). Elevated CO2 reduced the overall mineral concentrations 
in crops (-7.2%, -8.6 to -5.6); wild (-9.7%, -1 1 .6 to -7.8), herbaceous (-7.5%, -8.7 to -5.6), and 
woody (-9.6%, -12.1 to -7.6) plants; foliar (-9.2%, -10.8 to -7.6) and edible (-6.4%, -7.8 to -5.1) 
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Figure 5. The effect of CO2 in artificial enclosures. Change (%) in the mean concentration of chemical elements 
of plants grown in chambers, greenhouses, and other artificial enclosures under eC02 relative to those grown at 
ambient levels. Average ambient and eC02 levels across all the non-FACE studies are 365 ppm and 732 ppm 
respectively. Other details are in the legends for Figures 2 and 3. 
DOI: 10.7554/eLife.02245.007 



tissues. Including grains (-7.2%, -8.6 to -5.6). The cereal specific declines in grains are as follows: 
wheat (-7.6%, -9.3 to -5.9), rice (-7.2%, -1 1 .3 to -3.1), and barley (-6.9%, -10.5 to -3.2) {Figure 8). 
This is notable because wheat and rice alone provide over 40% of calories to humans. 

Discussion 

The analysis of all the data shows that eC02 shifts the plant ionome toward lower mineral content; the 
mean change across all the 25 measured minerals is (-8%, -9.1 to -6.9) {Figure 2). This shift, however, 
is hidden from low-powered statistical tests {Figure 1). Attaining adequate meta-analytic power 
reveals that the shift is: 

1. Empirically robust — evident in both artificial (chambers, greenhouses) and field (FACE) conditions 
{Figures 5 and 6). 

2. Geographically pervasive — ^found in temperate and subtropical/tropical regions {Figures 7 and 8). 

3. Systemic — affecting herbaceous and woody plants, crops, and wild plants, photosynthetic and ed- 
ible tissues, including wheat, rice, and barley grains {Figures 3, 4 and 8). 

Elevated CO2 alters plant C:N:P:S stokhiometry 

Not only does eCOa reduce the plant mineral content, but it also alters plant stoichiometry. Specifically, 
the effect of eC02 on N is nearly twice as large as its mean effect on minerals. The differential effect of 
eC02 on N (15%), and P (9%) and S (9%) translates into a -7% reduction in the plant N:P and N:S. 
In contrast to the lower N and mineral content, eC02 increased C content by 6% {Figures 2, 3 and 5). 
It follows then that eC02 increases C:P and C:S by 16%, and C:N by 25% confirming the previous findings 
of 19-27% higher C:N in plants grown in eC02 {Poorter et al., 1997; Stiling and Cornelissen, 2007; 
Robinson et al., 2012!). 
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Figure 6. The effect of CO2 at FACE centers. Change (%) in the mean concentration of chemical elements of plants 
grown in Free-Air Carbon dioxide Enrichments (FACE) centers relative to those grown at ambient levels. Average 
ambient and eC02 levels across all the FACE studies are 376 ppm and 560 ppm respectively. Other details are in 
the legends for Figures 2 and 3. 

DOI: 10.7554/eLife.02245.008 



Data scarcity 

The current dataset (available at Dryad depository) suffices to show the overall shift In the plant 
ionome. However, it would require much richer datasets to quantify differences among the shifts 
of various minerals and to assess shifts in the ionomes of individual species. Unfortunately, funding 
hurdles for analyzing fresh and archived samples harvested at FACE centers have significantly delayed 
progress In this area. Only two CO2 studies report selenium (Se) content {Hogy et a/., 2009, 2013), 
and none report data on tin (Sn), lithium (Li), and most other trace-elements. For many of the world's 
popular crops, pertinent data are non-existent or very limited, including (in the descending order of 
calories provided to the world's population, FAO, 2013): maize (the top C4 crop), soybeans (including 
oil), cassava, millet, beans, sweet potatoes, bananas, nuts, apples, yams, plantains, peas, grapes, 
rye, and oats. 

The current data scarcity, however, should not detract our attention from what is likely to be the 
overarching physiological driver behind the shift in the plant ionome — the C02-induced increase in 
carbohydrate production and the resulting dilution by carbohydrates. Let us take a closer look at this 
nutritionally important issue. 

TNC:protein and TNC:mmerals respond strongly to elevated CO2 

Carbohydrates in plants can be divided into two types: total structural carbohydrates (TSC; e.g., cellulose 
or fiber) that human body cannot digest, and total non-structural carbohydrates (TNC), most of 
which- — including starch and several sugars (fructose, glucose, sucrose, and maltose) — is readily 
digestible and absorbed in the human gut. Hence, for humans, TNC carries the most of caloric and 
metabolic load of carbohydrates. Out of the two types of carbohydrates, eCOa affects stronger the 
latter, boosting TNC concentration by 10-45% (Stiling and Cornelissen, 2007; Robinson et a/., 
2012). Furthermore, eC02 tends to lower protein in plant tissues {Taub et a/., 2008). Hence, we 
can reason that eC02 should exacerbate the inverse relationship found between TNC and protein 
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Figure 7. The effect of CO2 at various locations and latitudes. Locations of the FACE and Open Top Chamber 
(OTC) centers, which report concentrations of minerals in foliar or edible tissues, are shown as white dots inside 
colored circles. The area of a circle is proportional to the total number of observations (counting replicates) 
generated by the center. If the mean change is negative (decline in mineral content), the respective circle is blue; 
otherwise, it is red. The figure reflects data on 21 minerals in 57 plant species and cultivars. The shaded region 
(between 35 N and S latitudes) represents tropics and subtropics. 
DDI: 10.7554/eUfe.02245.009 



{Poorter and Villar, 1997). Considering thatJNC and protein are two out of the three primary macro- 
nutrients (with fats/lipids being the third), it becomes imperative to quantify changes in TNC:protein, 
when estimating the impact of altered plant quality on human nutrition in the rising CO2 world. 

Regrettably, TNC:protein is rarely reported by CO2 studies; instead C:N is used as a yardstick for 
accessing changes in the plant quality. However, C:N poorly correlates with TNCiprotein because protein 
is more C-rich than carbohydrates (C content in protein is 52-55% cf. 40-45% in carbohydrates). Thus, a 
h/gher carbohydratetprotein results in a lower C content. This means that C02-induced changes in 
nutritionally and metabolically important ratios — TNC:protein and TNC:minerals — can substantially 
exceed the respective changes in C:N. We can calculate changes in TNCiprotein using reported changes 
in TNC and protein (see 'Formula for calculating percentage changes in TNC:protein and TNCtminerals' 
in 'Materials and methods'). Table 1 compares C02-induced changes in C:N with respective changes in 
TNCiprotein. It shows that eC02 can elevate TNCiprotein up to fivefold higher than it does CiN. 

How shifts in TNCiprotein affect human nutrition is still unknown. New evidence, however, chal- 
lenges "the notion that a calorie is a calorie from a metabolic perspective" by showing that changes 
in dietary carbohydrateiproteinifat ratios affect metabolism and weight gain in humans {Ebbeling et 
a/., 2012). The new evidence supports an emerging view that while obesity is quantified as an imbal- 
ance between energy inputs and expenditures (Ha// et a/., 2011), it could also be a form of malnu- 
trition (We//s, 2013), where increased carbohydrateiprotein {Simpson and Raubenheimer, 200S) and 
excessive carbohydrate consumption {Taubes, 2013) could be possible culprits. 

Absolute CO2 effect on TNC. Spoonful of sugars for everyone? 

The baseline TNC content in plant tissues varies widely. In grains and tubers, it is very high, 50-85% 
of dry mass (DM). Therefore, in these tissues a percentage increase in TNC is arithmetically limited 
(e.g., a 60% increase is impossible). However, even a modest percentage increase in TNC-rich tissues 
can be nutritionally meaningful in absolute terms. For example, the FACE study of Porteaus et al. 
(2009) reports a 7-8% increase in starch concentrations in wheat grains, which translates to ~4 g 
of additional starch per 100 g DM. In contrast to grains and tubers, the baseline TNC level in photo- 
synthetic tissues is small (usually <25%), which makes large TNC increases possible. For example, 
Teng et al. (2006) reports that eC02 increased TNC by 76% in leaves of Arabidopsis thaliana. What is 
interesting here is that in absolute terms (per 100 g DM) the ~5 g TNC increase in Arabidopsis thaliana 
is comparable to the ~4 g TNC increase in wheat grains. 
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Figure 8. The systemic aspect of the CO2 effect. Change (%) in the mean concentration of minerals in plants grown 
in eC02 relative to those grown at ambient levels. All the results in the figure reflect the combined data for the 
foliar and the edible tissues. The number of total mean observations (m) for all the measured minerals across all 
the studies for each crop/plant group, experiment type, country, or region is shown with the respective statistical 
power Country specific and regional results reflect all the FACE and Open Top Chamber (OTC) studies carried in 
any given country/region. The number of total observations (with replicates) for all the minerals (not counting C and 
N) for each country is as follows: Australia (926), China (193), Finland (144), Germany (908), and USA (1156). Other 
details are in the legends for Figures 2 and 3. 
DOI: 10.7554/eLife.02245.010 



More generally, CO2 studies show that— Irrespective of the baseline TNC content — eCOa tends to 
boost TNC by a few grams (1-8 g) per 100 g DM of plant tissue {Poorter et al., 1997; Keutgen and 
Chen, 2001; Katny et al., 2005; Erfas et al., 2010; Azam et al., 2013). Note that such an infusion 
of carbohydrates into plant tissues, all else being equal, dilutes the content of other nutrients by 
-1-7.4%. Let us compare the dilution with its pragmatic and easily graspable analog — adding a 
spoonful of sugar-and-starch mixture. Table 2 shows that the CO2 effect on TNC:protein and 
TNC:minerals is stoichiometrically similar to the effect of adding a spoonful of carbohydrates to every 
1 00 g DM of plant tissue. 

Clearly, adding a spoonful of sugar sporadically to one's diet is not a cause for concern. However, 
the inescapable pervasiveness of globally rising atmospheric CO2 concentrations raises new questions: 
What are health consequences, if any, of diluting every 1 00 g DM of raw plant products with a spoonful 
of starch-and-sugar mixture? What are the consequences if the dilution is not sporadic but unavoidable 
and lifelong? These questions are better left for nutritionists, but it is worth noting that WHO (2014) 
conditionally recommends that intake of free sugars not exceed 5% of total energy, which is equivalent 
to 5-8 teaspoons of sugar for a typical 2000-3000 kcal/day diet. 

Below, I shift focus on a direct consequence of the C02-induced increase in carbohydrate production — - 
the mineral decline in plant tissues, and explore its potential effect on human nutrition. 

Plant minerals and 'hidden hunger' 

'Hidden hunger' — stems from poorly diversified plant-based diets meeting caloric but not nutritional 
needs. It is currently the world's most widespread nutritional disorder {Kennedy et al., 2003; 
Welch and Graham, 2005). It lowers the GDP of the most afflicted countries by 2-5% and is partly 
responsible for their Third World status {WHO, 2002; Stein, 2009). A paradoxical aspect of 'hidden 
hunger' is that the minuscule amount of minerals, which a human body requires, could be provided 
easily and inexpensively — at least in theory — to all people in need by fortifying foods with minerals. 
However, in practice, such required mineral levels do not reach large parts of the world's community. 
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Figure 9. Testing for publication bias. A funnel plot of the effect size (the natural log of the response ratio) plotted 
against the number of replicates/sample sizes (n) for each study and each mineral in the dataset for C3 plants. The 
plot provides a simple visual evaluation of the distribution of effect sizes. The blue line represents the mean effect 
size of eC02 on mineral concentrations: the decline of 8.39% (yielding the decline of 8.04% when back transferred 
from the log-form). The symmetrical funnel shape of the plot around the mean effect size indicates the publication 
bias in the dataset is insignificant {Egger et a/., 1997). 
DOI: 10.7554/eLife.02245.011 



The case of iodine Is illustrative: although iodized table salt nearly wiped out iodine deficiency in the 
industrialized world, a billion people still have no regular access to it, making iodine deficiency the 
leading cause of preventable brain damage, cretinism, and lower IQ in children {Welch and Graham, 
1999; WHO, 2002). Hence, the reality of logistic, economic, and cultural hurdles for fortification leaves 
the natural and bioavailable mineral content in food, and in plants in particular, to be the major, and 
sometimes the only, consistent mineral supply for a large part of mankind [White and Broadley, 2009; 
Bouis and Welch, 2010). This supply, unfortunately, is suboptimal for human nutrition with some 
of the consequences outlined below. 

Every third person in the world is at risk of inadequate Zn intake with its deficiency substantially 
contributing to stunting, compromised immunity, and child mortality {Brown et al., 2001; UNICEF, 
2009). Iron deficiency affects at least 2 billion people and is the leading cause of anemia that increases 
maternal mortality {WHO, 2002; UNICEF, 2009). Millions are Ca, Mg, and Se deficient [Stein, 2009; 
White and Broadley, 2009), including some population segments of developed countries [Rayman, 
2007; Khokhar et al., 2012). Ironically, a person can be obese and mineral undernourished — the so 
called 'hunger-obesity paradox' [Scheier, 200S), for example the many homeless in the US who rely on 
"cheap and energy-dense but low-nutrient" foods [Koh et al., 2012). With every third adult in the 
world being overweight or obese (Keats and Wiggins, 2014), WHO ranks both mineral undernutrition 
and obesity among the top 20 global health risks [WHO, 2002; Hill et al., 2003; Stein, 2009). While 
the role of mineral deficiency in obesity is still unclear, intriguing links have been found between the 
lower blood serum concentrations of Ca, Cr, Fe, Mg, Mn, Se, Zn, and increased body mass index (BMI), 
with most of the findings appearing in the last decade [Singh et al., 1998; Martin et al., 2006; 
Arnaud et al., 2007; Garcia et al., 2009; Payahoo et al., 2013; Yerlikaya et al., 2013). 
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Roumet et al. (1999)* 
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Porteaus et al. (2009) 


wheat grain (high N) 
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Poorter et al. (1997) 


meta-analysis 
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Robinson et al. (2012) 


meta-analysis 


27 


39 
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C02-induced changes (%) in C:N (a quality indicator often used in CO2 studies) and in TNC:protein (a rarely used 
but nutritionally important indicator) for wheat grains and for foliar tissues of various plants. The results shows that 
in the same plant tissue, eC02 can increase TNC:protein up to several-fold > C:N. Significant C02-induced shifts in 
the ratio of major macronutrients are probable. Hence, it is important for CO2 studies to start accessing and 
reporting changes in TNCiprotein. 
*in lieu of protein, N content is used. 
DOI: 10.7554/eLife.02245.012 



How can the C02-induced depletion of minerals in crops affect humans? I emphasize that the 
impact of C02-induced shifts in the quality of crops on human health is far from settled. The purpose 
of what follows is not to make definitive claims but to stimulate research into this important but 
unresolved issue. 



Table 2. Comparing the effect of CO2 to the 
effect of adding 'a spoonful of sugars.' 

Effect of 



Plant quality 
indicator 


Effect of adding 
5g of TNC (%) 


elevated 
CO2 (%) 


Grains and tubers: 


TNC 


2.6 


1 to 15 


TNC:protein 


7 


6 to 47 


TNCiminerals 


7 


6 to 28 


protein 


-4.8 


-14to -9 


minerals 


-4.8 


-lOto -5 


Foliar tissues: 


TNC 


27 


15 to 75 


TNC:protein 


33 


26 to 125 


TNC:minerals 


33 


24 to 98 


protein 


-4.8 


-19to -14 


minerals 


-4.8 


-12 to -5 



Changes (%) in various plant quality indicators caused by: 
(1) Adding a teaspoon of TNC (~5g of starch-and-sugars 
mixture) per 1 0Og of dry mass (DM) of plant tissue, 
an:d (2) growing plants in twice-ambient CO2 atmosphere. 
Changes due to the addition of TNC are calculated 
assuming:the baseline TNC content of 65% for grains 
and tubers, and 1 5% for foliar tissues. The C content is 
assumed to be -42% for plant tissues and TNC. 
DOI: 10.7554/eLife.02245.013 



Stokhiometrk thought experiment 

A randomized controlled trial for a human diet 
based exclusively (directly or indirectly) on plants 
grown in eC02 is unlikely and ethically question- 
able; and even if feasible, the trial might take 
years to generate results. In lieu of relevant data, 
we can employ a thought experiment. While such 
'experiments' are usually reserved for physical 
sciences, any living system, notwithstanding its 
complexity, adheres to simple but irrefutable 
elemental mass balance, which can help us to 
elucidate plausible scenarios. 

For simplicity, let us focus on one question: 
how can a 5% reduction in the plant mineral 
content affect human nutrition? Thus, we ignore 
other potential or likely CO2 effects: for example 
higher agricultural yields; altered concentrations 
of lipids, vitamins, and polyphenols; substantially 
higher TNC:protein and TNC:minerals; differential 
C3 and C4 plant responses; changes in the phytate 
content that affects mineral bioavailability {Manoj- 
Kumar, 2011); and multiplicative health effects of 
the concomitant declines of many minerals in the 
same tissue. 

Suppose that starting tomorrow and without 
our knowledge, the baseline mineral content of 
all plants on Earth drops by 5%. A self-evident 
but easily overlooked mass-balance law tells us 
that neither thermal nor mechanical processing of 
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raw plants enriches them with minerals (i.e., transmutations are impossible). Thus, the mineral decline 
in raw crops will follow into plant-based foods (except for a few food items that are fortified with 
certain minerals in some countries). 

We can safely assume that the individuals, whose dietary intake of each essential mineral has 
exceeded the recommended dietary intake (RDI) by >5%, will be unaffected by the depletion. This 
leaves us with the majority of the human population, whose diet is either at risk of deficiency or already 
deficient in atleast one mineral (WHO, 2002; Kennedy et al., 2003; Stein, 2009). Though a human 
body can synthesize complex compounds (e.g., vitamins K and D, non-essential amino acids), the mass 
balance low implies that no organism can synthesize any amount of any mineral. Therefore, to com- 
pensate for the mineral deficit, an organism has to increase mineral intake (or, otherwise, endure the 
consequences of the deficit). Taking supplements or intentionally shifting one's diet toward mineral- 
rich foods, for example animal products, can eliminate the deficit. Such dietary changes, however, 
presuppose behavioral adjustments on the part of the individuals who are aware of their mineral deficiency 
and have both the means and motivation to address it. A simpler way to compensate for the mineral 
deficit is to increase food intake, whether consciously or not. (The notion of compensatory feeding is 
not entirely alien — herbivores do increase consumption by 1 4-1 6%, when consuming plants grown 
in eC02; Stiling and Cornelissen, 2007; Robinson et al., 2012). 

For a calorie deficient person, eating 5% more (to be exact 5.26%, because 1.0526*. 95 » 1) is likely 
to be beneficial. However, for a calorie sufficient but mineral deficient person, eating 5% more could 
be detrimental. The dynamic mathematical model of human metabolism, which links weight changes 
to dietary and behavioral changes {Hall et al., 201 1), can help to quantify the effect of a prolonged 
5% increase in food intake. When parameterized with anthropometric data for an average moderately 
active American female (age 38, height 163 cm, weight 76 kg, BMI 28.6, energy intake 2431 kcal/day 
[1 01 71 kJj) {Fryar et al., 2012; CIA, 2013), the model outputs a weight gain of 4.8 kg over a 3-year period, 
provided all other aspects of behavior and diet remain unchanged. For a male, the respective weight 
gain is 5.8 kg. The results are congruent with Hill et al. (2003), who argued that a 4-5% difference in total 
daily energy intake, a mere 100 kcal/day, could be responsible for most weight gain in the population. 

The above 'experiment' suggests that a systemic and sustained 5% mineral depletion in plants can 
be nutritionally significant. While the rise in the atmospheric CO2 concentration is expected to be 
nearly uniform around the globe, its impact on crop quality might unequally affect the human population: 
from no detrimental effects for the well-nourished to potential weight gain for the calorie-sufficient 
but mineral-undernourished. 

Has rising CO2 already altered the plant ionome? 

The rise in CO2 levels over the last 18-30 years has already been implicated in the two effects that can 
influence the plant ionome: higher C assimilation and plant growth {Donohue et al., 2013), and lower 
transpiration {Keenan et al., 2013). Considering that over the last 250 years, the atmospheric CO2 concen- 
tration has increased by 1 20 ppm — an increase that is not far from the mean 1 84 ppm enrichment in the 
FACE studies — it is plausible that plant quality has changed. Indeed, declines in mineral concentrations 
have been found in wild plants and in crop fruits, vegetables, and grains over 22-250 years {Penuelas and 
Matamala, 1993; Duquesnay et al., 2000; Davis et al., 2004; Ekholm et al., 2007; Fan et al., 2008; 
Jonard et al., 2009). While the mineral declines in crops can be an unintended consequence of the 
Green Revolution that produced high-yield cultivars with altered mineral content {Davis et al., 2004; 
Fan et al., 2008), the reason for the mineral declines in wild plants cannot be attributed to it. 

Can eC02 directly affect human health? Hersoug et al. (2012) proposed that rising CO2 promotes 
weight gains and obesity in the human population directly (via breathing) by reducing the pH of blood 
and, consequently, increasing appetite and energy intake. Weight gain has been observed in wild 
mammals, lab animals, and humans over the last several decades {Klimentidis et al., 2011). However, 
it is not clear what role, if any, the rising CO2 could have played either directly (breathing) or indirectly 
(altered plant quality). And disentangling the rising CO2 effect from other plausible factors currently 
does not seem feasible due to scarce data. This brings us to the broader issue of detecting — amid high 
local noise — signals that are small in their magnitude but global in their scope. 

Hidden shifts of global change 

While some scientific areas (e.g., genomics, bioinformatics) have experienced a data deluge, many 
areas of global change, including the issue of shifting plant quality, have been hindered by chronic 
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data scarcity. Fortunately, researchers worldwide have been steadily generating data on the effects of 
eCOa on the chemical composition of plants. It is their collective efforts that have made it possible to 
reveal the C02-induced shift in the plant ionome. 

Human activities profoundly alter the biogeochemical cycle not only of C but also of N, P, and S, 
which are central to all known life forms. It is plausible that other subtle global shifts in the physiology 
and functioning of organisms lurk amid highly noisy data. The small magnitude of such shifts 
makes them hard to detect and easy to dismiss. But by virtue of being global and sustained, the shifts 
can be biologically potent. Revealing hidden shifts requires plentiful data to attain sufficient statistical 
power. (For example, Rohde et a/. (2013) analyzed 1 4 million mean monthly local temperature records 
to uncover the 1.5°C rise in the global average temperature since 1753 — undoubtedly a potent but a 
very small change relative to the variations of tens of degrees in local temperature.) 

New data on the effects of eCOa on plant quality (e.g., minerals, TNC: protein, TNCrminerals, lipids, 
bioavailability of nutrients) can be generated very cost-efficiently by analyzing fresh and archived plant 
samples collected at FACE centers worldwide (the project leaders of many centers are keen to share 
such samples; PS Curtis, BA Kimball, R Oren, PB Reich, C Stokes; IL personal communication, July, 
2006). With regard to minerals, the application of the high-throughput techniques of ionomics {Salt 
et al., 2008) can generate rich phenotypic data that can be linked with functional genomics. Such 
analyses will shed more light on changes in plant quality in the rising CO2 world. Anticipating and 
assessing such changes will help not only in mitigating their effects but also in steering efforts to breed 
nutritionally richer crops for the improvement of human health worldwide. 

Materials and methods 
Search for data 

I searched Google Scholar, Google, PubMed, the IS! Web of Science, AGRICOLA, and Scopus to find 
relevant articles with sensible combinations of two or more of the following search-words: elevated, 
rising, CO2, carbon dioxide, ppm, FACE, effects, content, concentration, %, mg, dry matter, micronu- 
trients, plant(s), crop(s), tree(s), C3, C4, foliar, leaves, grains, seeds, tubers, fruits, minerals, chemical 
elements, and names/symbols of various chemical elements (e.g., zinc/Zn). I found additional studies 
from references in the articles identified in the initial searches. 

Study suitability and data selection criteria 

Among all plant tissues for which mineral concentrations are reported in the literature, the most abun- 
dant data are on foliar tissues (leaves, needles, shoots), and — for herbaceous plants — on above ground 
parts. Hence, focusing on the foliar tissues and above ground parts allows one to maximize the number of 
independent observations of the effect of eCOa on each mineral. Although the data on edible parts of 
crops are scarcer, a dataset on crop edible tissues was compiled due to their direct relevance for 
human nutrition. 

The following objective and uniform criteria were applied for deciding which studies to include into 
the dataset: (1) a study grew plants at two or more CO2 levels, (2) a study directly measured the con- 
tent of one or more minerals in foliar or edible plant tissues at low (ambient) and high (elevated) CO2 
levels, and (3) a study reported either absolute concentrations at each treatment or relative change/ 
lack thereof in the concentrations for each mineral between treatments. Studies that indirectly deduced 
mineral concentrations, reported data on N but not on any mineral, exposed only a part (e.g., a branch) 
of the plant, used super-elevated or uncontrolled CO2 levels were not included. Table 3 lists all the 
studies together with their respective species/cultivars and CO2 enrichment levels (the dataset with all 
the details is deposited at Dryad and GitHub). When a study reported the low CO2 level as 'am- 
bient' with no specific numerical values, then I used the Keeling curve to approximate the ambient 
CO2 level for the year the study was carried out. 

The following data-inclusion rules were applied to the studies with multiple co-dependent datasets 
for the foliar dataset: (1 ) the lowest and the highest CO2 levels for studies with multiple CO2 levels, 
(2) the control and single-factor CO2 for studies with environmental co-factors (e.g., observations from 
combined eC02 and ozone experiments were excluded), (3) the highest nutrient regime when the 
control could not be identified in a study with multiple nutrient co-factors, (4) the last point, that is the 
longest exposure to ambient/eC02 for studies with time series, (5) the most mature needles/leaves 
for studies reporting foliar tissues of various ages. If, in rare instances, a publication reported three 
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effects on the plant 


1 ll^fl 1 IC 




Common 














t*rop 


+C02 


Country 


IXt. It.1 CIIWC 




r 1 1 dUlc Lice 


No 


260 




Ii^vel uieCK/ f rrO 


ZA/^^f ft ir^ft itrt^ 

/-\cfc?f luL/iuin 


rcU rlldUifc; Lice 


No 


200 


USA 






grass 


No 


340 


UK 








No 


250 




IVcWucf cr alt, 1 7 7D 


/-\iIIUb y/L/L/ncjoa 


cUvjcl Lr cc 


No 


350 


UK 


T^mnaw+nn at of 9/1/1? 
ISIJipSltOn CL altg £U\IO 


/-Kiljl il LUi ila LJfcrL/lc;/ 


rdlliiUicbL Lice 


No 


440 




IxanQWanf, ^KiU 1 


/-\t 1 iDS ijhta (JIUIJiDoa 




No 


180 


USA 


U/Micmon a+ of 9/119 
nousijlan CI alt, £\I l£ 


AfQlDidopsis thsllQIlQ 


LI Id Ic CI ebb 


No 


450 




Mill a+ of 9/11 '7 
IVfU eZ 31. f £.\1 \J 


aLJiLJL/LJbli} Liiallaiia 


LI Id Ic CI ebb 


No 


330 






DcLUIa L/ciUJUla i\DL\\ 


1 ri + res ci 
LJ 1 1 Cl 1 LI cc 


No 


349 


Finland 


^freonan at of 9/1/m 




Qrass 


No 


230 




P/^ffai/ at of 9/111 


Dl UI 1 lUb LfcrCt-QfUffJ 


CI led Ly r dbb 


No 


150 




RIanlf at of 9/lflA 
DlaflK SL alt, £,\J\JO 


Dl UI nUb LfcrCtL'fUlM 


CI led Ly 1 dbb 


No 


150 




Rfonir at of 9/111 
DId/IK er alt, £\J 1 1 


K^aiiLii la VUlLJal lb 


lledLllcl blilLJL) 


No 


200 




Wnrtrltn at of 1009 
VVOOOin ex ait, l Tr£ 


I ^ * c ^ 3 *o 3 n c f c 
t^ciL/o Ldf (aUfcrtlb/b 


r eCl UUCj LI cc 


No 


200 


USA 


pMwfr! at of Of\n^ 

nnzi St alt, 1 


11 ybai tumniUlii 1 1 liji iiUIIUI 1 1 


CI II ybd 1 1 LI 1 


No 


325 




1^1 lafimr at of 1001 

iKUenny ei ait, 1 77 1 


K^Ul 1 lUb llUilUa 


Ciuyvvc/uci Lice 


No 


200 


USA 


PSwfri at of 9/1/11 

nnzi et ait, £iju i 


laLJUb byiVaLK^a 


LJccCil Lice 


No 


260 




^i/at-^ia/*l^ 1 0O? 


laLJUb byiVaLli^a 


LJccCIl Lice 


No 


300 




l?/^Wanl^it*/'fian at of 9/1/10 

nooenKircnsn st ait, 


icbLUl^a L)l aLci ibib 


iiicdciuvv icbCLic 


No 


320 




^veioiecK, 1 77*3 


ntratUCd ViViL/al a 




No 


340 


UK 


l3ovtar at of 
oaXtSr alt, 1 774 


rnfiUc/ofa L/laytcyalia 


laiiiiL/icoL Lice 


No 


440 




IxanoWanf, ^w\J 1 




P 1 1 1 c rv^ 1 ^ 
CIIIULL b IIIIIKLJcd 


No 


325 


USA 


Lfi.MMota at of OflFiA. 

nunyaLe ex ait, £\iu*t 


i—dl IK KdcrfilUifcrf / 


IdlCII Llec 


No 


335 


Japan 


Cf^fnono at of 9/1/17 

oninano ex ait, £.u\i i 


I—CL/ILIIUI 1 1 IciLIIUIIUni 


peppcl y 1 dbb 


No 


339 




DianK ana uerner, ^uuh 


i—IUUIUan lljal oLV/dC/nUd 


b Wee Ly LI 1 T 1 Lice 


No 


200 


USA 


PSwfri at of 9/1/11 

nnzi ex ait, £iju i 


i—IUUIlJal 1 ILJai biyi al^UtUcI 


SVVCCLyLJIII Lice 


No 


167 


USA 


Inhncnn at of 9/1/1^ 

tjonnson ex ait, £\iuh 


i—ILiUILJaniLJal b tyf dCi f i Ud 


b Wee Ly LJ 1 11 Lice 


No 


1 56-200 


USA 


fVfotofi at of 9/1/10 


l_Jf ((JUfcrf iCJi Qf 1 LUIlLJiIcI a 


LLI III-' LI ce 


No 


325 




O'fVfaiff at of ^0fl7 
K,/ iveill ex alt, IrO/ 


I—UIIUIII LJcIclllic 


grass 


No 


320 




KjveraiecK., i rro 


LUMU/fl Uclclillc 




No 


290 


Germany 


C/*lianf^ at of 1007 


i—UL)ll lUb aiL/Ub 


\Air~\ |■^d 1 ir^ 1 n 

Will LC 1 LJLJ 1 1 1 


No 


550 




^ommf^aff on^ Co^a 9/1/19 
\^alTipDeil ano Oaye, £\J\I£ 


Lycium pailidum 


shrub 


No 


180 


USA 


Housman et a/., 2012 


Nephrolepis exaitata 


fern 


No 


650 




Nowak et a/., 2002 


Pelargonium x hortorum 


geranium 


No 


330 




M/shra eta/., 2071 


IVIdVtiriCK VVlllLti 












r/Ccrd aUicb ixdibL. 


bUI UCc LI cc 


No 


350 




rTirrmann et ai*, i yyo 


r/Ccrd dU/co IxdibL. 


bpi UCc LI tic 


No 


300 




nouennifcnen er af</ ^i/U7 


r/Ccrd dLJ/cb IxdibL. 


bprucc Ulcc 


No 


300 






r/Ccrd fULJcrflb 


bUlLJCc Lice 


No 


350 




Cfiinf^i/ a+ of 1009 

onipiGy ex ait, 1 77^ 


Pir\iic f~tf~tr\/~lOKi~iC^ 
1 IIILlJ LJt-/l ILJCi LJja 


i— J IJ 1 C Lice 


No 


346 


USA 


WalKcf cl all/ 


Pinus ponderosa 'Laws.' 


pine tree 


No 


350 


USA 


Johnson et a/., 1997 


P/nus sy/vestr/s 


pine tree 


No 


331 




Luoma/a et a/., 2005 



Table 3. Continued on next page 
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Table 3. Continued 



Common 



SpGCIGS 


ndlTIG 


Crop 




i^ouniry 


ReTG re n ce 


Pinus syivGstris 


pine tree 


INO 


ZZD 


Finland 


utnainen et ai., £.\jUU 


Pinus tBQcis 


loblolly pine tree 


INO 


9nn 
zuu 


UoA 


rinzi et ai.f ^uu 1 


Pinus taec/a 


pine tree 


INO 


9nn 
zuu 




niatau et ai., zuuy 


Poa 3!pin3 


grass 


INO 


O^U 


UIn 


tsaxTer et ai.f 1774 


PoB dlpitld 


grass 


INO 


OM-U 


U In 


DdXTer et ai.f iw/ 


PtGndiurn dcjuilinuni 


fern 


INO 


OZU 




^neng et ai., ^uuts 


Ptsr/c//um rGvoiutum 


fern 


INO 


OZU 




^neng et ai./ ^uuts 


Pteris vitt3t3 


fern 


INO 


oZU 




^neng et ai., zuuo 


CluQfcus chdpmBnii 


oak tree 


INO 


oOU 




fVdtdff et ai./ ^L/uy 


QuGTcus gemindtd 


oak tree 


INO 


jOU 


UoA 


jonnson et ai., iduuo 


QuGTCUS QQtDindtS 


oak tree 


INO 


jOU 


1 KA 
UoA 


fvataff et ai., ici/uy 


QuGrcus rnyrtifoUd 


oak tree 


INO 


jOU 


I IQA 
UoA 


jonnson et ai., <cuuo 


QuGrcus myrtifolia 


oak tree 


No 


350 


USA 


Natali et a/., 2009 


QuGrcus suber 


cork oak tree 


No 


350 




Niinemets et ah, 1999 


Schizachyrium scopahum 


grass 


No 


230 




Polleyet a/., 2011 


Sorghastrum nutans 


grass 


No 


230 




PoUeyet a/., 2011 


Sporobolus kGntrophyllus 


grass 


No 


330 




Wilseyetai, 1994 


Thfoiium alGxandhnum 
'Pusa Jayant' 


berseem clover 


No 


250 


India 


Pal et a/., 2004 


Trifolium pratGnsG 


red clover 


INO 


oZU 




(.^veraiecK/ lyyj 


Trifoliutn repens 


white clover 


INO 


oZU 




cverafecK, / yyj 


TrifoliuiD repens 


white clover 


INO 


9on 
zyu 


Germany 


C^UrtnL- n4i -.1 •1007 

^cnenfc et ai., lyy/ 


TrifoliuiD fGpGns 


white clover 


INO 


0 1 D 




1 fan et ai 1 4 


Trifolium repens 'Regal' 


white clover 


INO 


ooU 




neagie et ai., iwo 


ValiisnGria spinulosa 


macrophyte 


No 


610 




Yan et a/., 2006 


Apium gravGoiGns 


celery 


Yes 


670 




Tremblay et ai, 1988 


Brassica juncGa 'Czern' 


mustard 


Yes 


500 


India 


Singh et al., 2013 


Brassica napus 'Qinyou 8' 


rapeseed 


Yes 


615 




Tian et al., 2014 


Brassica napus 'Rongyou 10' 


rapeseed 


Yes 


615 




Tian et al., 2014 


Brassica napus 
'Zhongyouza 12' 


rapeseed 


Yes 


615 




Tian et al., 2014 


Brassica napus Campino 


oilseed rape 


les 


1 OA 
1 UO 


Germany 


nogy et ai., £U i u 


Brassica rapa Grabe 


turnip 


les 


Ann 
ouu 




Azani et ai., £U u 


Citrus aurantium 


orange tree 


les 


"^nn 

oUU 


I IQA 
UoA 


Penuelas et al., 1997 


Citrus madurGnsis 


citrus tree 


les 


Ann 
ouu 




fveutc^en dfiu K^nen, ^cuu i 


Cucumis sativus 


cucumber 


Yes 


650 




Peet etal., 1986 


Daucus carota 'T-1-1 1 1 ' 


carrot 


Yes 


600 




Azam et al., 2013 


Fragaria x ananassa 


strawberry 


Yes 


600 




Keutgen etal., 1997 


GlycinG max 'Mem' 


soybean 


Yes 


360 


USA 


Prior et al., 2008 


GiycinG max 'Mem' 


soybean 


Yes 


200 




Rodriguez et al., 201 1 


Gossypium hirsutum 
'Deltapine 77' 


cotton 


Yes 


180 


USA 


Huluka et al., 1994 


HordGum vuigarG 


barley 


Yes 


175 


Germany 


Erbs et al., 2010 


HordGum vulgarG 'Alexis' 


barley 


Yes 


334 


Germany 


Manderscheid et al., 7995 



Table 3. Continued on next page 
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Table 3. Continued 



Common 



Sp6Ci6S 


ndlTIG 


Crop 




^ounxry 


ReTG re n ce 


Hordeurn vuigBrs 'Arena' 


barley 


les 




Germany 


iviancierscneici et ai,, iwd 


HordQurn vuiQdfG Europa 


barley 


les 


Ann 




U44f A 4f on/is 


Hordeurn vufgdre 'Iranis' 


barley 


les 






Perez-Lopez et al., 2014 


Hordeum vu/^ars Therssa 


barley 


les 


1 70 

1 /u 


Germany 


wrODiewnz et ai., z*} 


Lsctucs S3tiv3 BRM 


lettuce 


les 


oUo 






Ldctucd sdtivd Mantilla 


lettuce 


les 








Lactuca sativa 'MV 


lettuce 


Yes 


308 




Baslam et a/., 2012 


Lactuca sativa 'Waldmann's Green' 


lettuce 


Yes 


600 




McKeehen et ai, 1996 


Lycopersicon escuientum 'Astra' 


tomato 


Yes 


600 




Khan et a/., 2013 


Lycopersicon escuientum 'Eureka' 


tomato 


Yes 


600 




Khan et a/., 2013 


Lycopersicon escuientum 'Mill.' 


tomato 


Yes 


360 




L/ et a/., 2007 


Lycopersicon escuientum 
/.neza ouv 


tomato 


Yes 


450 




J/n et a/., 2009 


Mangifera indica Kensington 


mango tree 


les 






dcnaTTer ancf wniiey, lyy/ 


Mangifera indica 'Tommy Atkins' 


mango tree 


les 






>cfiaTTer ancf wniieyf lyy/ 


Medicago sativa 


airaiTa 


les 


0 1 D 




1 fan et ai., 14 


Medicago sativa Victor 


dl 1 dlTd 




1 uu 


1 Ik* 

U IN 


Af-nawany et ai., ^cu lo 


Oryza sativa 


rice 


les 


zuu 


China 


rang et ai., <cuud 


Oryza sativa Akitakomachi 


rice 


les 


ZU3 — zou 


Japan 


LieTTenng et ai., ^u(/4 


Oryza sativa 'Akitakomachi' 


rice 


les 




Japan 


Y^makawa et a/., 2004 


\jryza saxwa DKrxiunan ov 


rice 


les 


91 n 
Z 1 u 


Bangladesh 


Kazzaclue et ai., ^L/uy 


Oryza sativa 'Gui Nnong Zhan' 


rice 


les 


ouu 




LI et ai.f Zulu 


Oryza sativa IR 72 


rice 


les 


90A 

zyo 


— — 

Philippines 


7iel/-a ii* al -1007 

^isKa et ai., i yy/ 


Oryza sativa 'Japonica' 


rice 


les 


9nn 
zuu 


China 


jia et ai., tcuu/ 


Oryza sativa Jarrah 


rice 


les 


oDU 




Seneweera and Conroy, 1997 


Oryza sativa Khaskani 


rice 


les 


91 n 
Z 1 u 


Bangladesh 


Kazzacfue et ai.f ^uuv 


Oryza sativa 'Rong You 398' 


rice 


les 


jUU 




Li et ai., icuiu 


Oryza sativa Shakkorkhora 


rice 


les 


91 n 

Z 1 u 


Bangladesh 


Kazzaclue et ai.f ^uuv 


Oryza sativa Shan You 428 


rice 




ouu 




1 S £X4- StI 9/11/1 

LI et ai.f lu 


Oryza sativa Tlan You 390 


rice 


les 


ouu 




1 I al 9/11/1 

LI et ai., f u 


Oryza sativa Wu Xiang jing 


rice 


les 


9nn 
zuu 


China 


ulio et ai., £U 1 1 


Oryza sativa 'Wuxiangjing 14' 


rice 


les 


9nn 
zuu 


China 


Aili -.1 9/1/17 

fwa et ai., ^uu/ 


Oryza sativa Wuxiangjing 14 


rice 


les 


9nn 
zuu 


China 


Vanji ^+ al 9/1/17 

Tang et ai., £.\ju/ 


Oryza sativa 'Yin Jing Ruan Zhan' 


rice 


les 


ouu 




Li et ai.f Zulu 


On/za sativa 'Yue Za 889' 


rice 


Yes 


500 




Liet al.. 2010 


Phaseolus vulgaris 'Contender' 


bean 


Yes 


340 




Mjwara et al., 1996 


Phaseolus vulgaris 'Seafarer' 


bean 


Yes 


870 




Porter and Grodzinski, 1984 


Raphanus sativus 'Mino' 


radish 


Yes 


600 




Azam et al., 2013 


Raphanus sativus 'Cherry Belle' 


radish 


Yes 


380 




Barnes and Mirrmann, 1992 


Raphanus sativus 'Giant White 
Globe' 


radish 


Yes 


600 




McKeehen et al., 1996 


Rumex patientia x R. 
Tianschanicus 'Rumex K-V 


buckwheat 


Yes 


615 




Tian et al., 2014 



Table 3. Continued on next page 
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Table 3. Continued 



Common 



Species 


name 


Crop 




i^ounxry 


ReTe re n ce 


Secale cereals 'Wintergrazer-70' 


rye 


les 


0 1 D 




1 fan et ai,, 1 4 


Soianum lycopersicum '76R MYC+'tomato 


les 


^on 

JYU 




K^avaQltBiO ex 31, f ^UU/ 


Soianum lycopersicum 'rmc' 


tomato 


Yes 


590 




Cavagnaro et a/., 2007 


Soianum tuberosum 


potato 


Yes 


500 




Cao and Tibbitts, 1997 


Soianum tuberosum 'Bintje' 


potato 


Yes 


170 


Germany 


Hogy and Fangmeier, 2009 


Soianum tuberosum 'Bintje' 


potato 


Yes 


278-281 


Sweden 


Piikki et a/., 2007 


Soianum tuberosum 'Bintje' 


potato 


Yes 


305-320 


Europe 


Fangmeier et a/., 2002 


Soianum tuberosum 'Dark 
Red Norland' 


potato 


Yes 


345 


USA 


Heagle et al, 2003 


Soianum tuberosum 'Superior' 


potato 








neagie et ai,f <cUUo 


Sorghum bicolor 


sorghum 


les 


oou 




rrlOr et ct/., £UUo 


Spinacia oleracea 


spinach 


les 




India 


jam et ai,, £,uu/ 


Trigonella foenum-graecum 


fenugreek 


les 




India 


jam et ai., icuu/ 


Triticum aestivum 


wheat 


les 


1 / D 


Germany 


ciDS et ai., £l/ i u 


Triticum aestivum 'Ningmai 9' 


wheat 


les 


zuu 


China 


fwa et ai./ zuu/ 


Triticum aestivum 'Triso' 


wheat 


les 


1 

\ 3U 


Germany 


nogy et a/., tcuuv 


Triticum aestivum 'Triso' 


wheat 


les 


1 3U 


Germany 


LJi;— „ -,1 9/1-1 
nogy et a/., lo 


Triticum aestivum 'Alcazar' 


wheat 


les 






cfe fa fuente et af., zuuu 


Triticum aestivum 'Batis' 


wheat 


les 


1 / u 


Germany 


wrODiewiTz et af., / j 


Triticum aestivum 'Dragon' 


wheat 


les 


one oofi 


Sweden 


P/e/je/ and Daniclssorig 2009 


Triticum aestivum 'HD-2285' 


wheat 


les 




India 


faf et af., ^uuo 


Triticum aestivum 'Janz' 


wheat 




1 oo 


Australia 


rernanao et af., zui4 


Triticum aestivum 'Jinnong 4' 


wheat 


les 


A1 

0 1 D 




I lan et af 1 4 


Triticum aestivum 'Minaret' 


wheat 


les 


97 J? 
Z/ O 


Germany 


FangmeiGr et a/., 1997 


Triticum aestivum 'Minaret' 


wheat 


les 


OUU 


Europe 


Fangmeier et a/., 1999 


Triticum aestivum 'Rinconada' 


wheat 


les 






ae fa fuente et af., i£uuu 


Triticum aestivum 'Star' 


wheat 


les 




Germany 


iviancierscneici et af., ivvd 


Triticum aestivum 'Turbo' 


wheat 


les 




Germany 


ivianoerscneiQ et af., fyys 


Triticum aestivum 'Turbo' 


wheat 


les 


jOU 




wu et af., ^i/u4 


Triticum aestivum 'Veery 10' 


wheat 




H 1 u 




^acfiefA af 9/1*19 

carifSfe et af., 1^ 


Triticum aestivum 'Yangmai' 


wheat 


Yes 


200 


China 


Guo et a/., 2017 


Triticum aestivum 'Yitpi' 


wheat 


Yes 


166 


Australia 


Fernando et a/., 2012a 


Triticum aestivum 'Yitpi' 


wheat 


Yes 


166 


Australia 


Fernando et a/., 2012b 


Triticum aestivum 'Yitpi' 


wheat 


Yes 


166 


Australia 


Fernando et a/., 2012c 


Triticum aestivum 'Yitpi' 


wheat 


Yes 


166 


Australia 


Fernando et a/., 2014 



The table provides species name, common name, the type of experimental set up, the level of CO2 enrichment, 
and indicates whether the species is a crop. Countries are listed only for FACE and OTC type experiments with 
'Europe' accounting for combined data from Belgium, Denmark, Finland, Germany, Sweden, and the UK. 
DOI: 10.7554/eLife.02245.014 



or more separate datasets for the same species or cultivar, the data were averaged prior to the inclu- 
sion into the foliar dataset. For the edible tissue dataset, the study inclusion rules were the same as for 
the foliar dataset with the following exception: due to relative scarcity of data for edible tissues, the 
data with co-factors were included in the dataset (e.g., observations from combined eC02 and ozone 
experiments were included). The 'Additional info' column in the dataset specifies exactly what datasets 
were extracted from each study with multiple datasets. 
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The above publication-inclusion and data-inclusion rules allow treating each study as independent 
in the dataset. At no instance, potentially co-dependent observations (e.g., multiple observations of 
the same plant throughout a growing season or observations of various parts of the same plant) were 
included in either the foliar or the edible dataset as separate studies. I used GraphClick v.3.0 and 
PixelStick v.2.5 to digitize data presented in a graphical form, for example bar charts. 

The foliar dataset covers 4733 observations of 25 chemical elements in 110 species and cultivars. 
The edible tissues dataset covers 3028 observations of 23 elements in 41 species and cultivars. The 
FACE studies cover 2264 observations of 24 elements in 25 species and cultivars. The two datasets 
reflect data on 125 C3 and 5 C4 species/cultivars. 

Effect size measure 

While the amount of statistical details provided in each study varies considerably, the following data were 
extractable from each study: (1) the relative change (or lack thereof) in the mean concentration between 
the low and the high CO2 treatments: {E-A)/A, where A and £ are the mean concentrations of an element 
at the low and the high CO2 treatments respectively, (2) the sample size or the number of replicates (n). 

Since a decrease in the concentration of a mineral is limited to 100%, but an increase in its 
concentration is theoretically unlimited, a standard technique was applied to reduce biases 
towards increases. Specifically, the natural log of the response ratio, that is ln(E/A), was used as 
the effect size metric (e.g.. Hedges et al., 1999, Jablonski et al., 2002; Taub et al., 2008). The 
response ratio, r = E/A, was calculated from the relative change as follows: r = 1+(E-A)/A. After 
performing statistical analyses, I converted all the results back from the log form to report them 
as ordinary percent changes. 

Making results replicable 

Published meta-analytic and biostatistical results need to be replicable and reproducible, and the pro- 
cess of replication needs to be made as easy as possible and clearly traceable to the original sources 
(Peng, 2009). In this regard, I have made the following efforts to ease the replication (from the original 
sources) of each and every result presented here: 

1. While copyright restrictions do not permit posting the original published data sources online, I 
will share, upon request, all the data sources in PDF form, where all the pertinent data are clearly 
marked for easy identification, thus removing any potential ambiguity about what data were ex- 
tracted from each study. 

2. The entire dataset for the foliar and the edible tissues is available at Dryad digital depository, www. 
datadryad.org, under 10.5061 /dryad. 6356f. The dataset is available as an Excel file (formatted for 
easy viewing) and as a 'CSV file; the latter is made-ready (tidy) for analysis with open-source (R 
Core Team, 2014) and commercial statistical packages (e.g., SPSS). 

3. An executable R code to generate individual results is available with the dataset at the above- 
mentioned depository and at GitHub: https://github.com/loladze/co2. Assistance for replicating 
any result and figure presented in this study will be provided to any interested party. 

Statistical analysis 

I performed all the analyses using R (R Core Team, 2014), SPSS v. 21 (IBM, Armonk, NY, USA) and 
G*Power 3 (Fau/ et al., 2007). Meta-analytic studies often weight effect sizes by the reciprocal of their 
variance, which tends to give a greater weight to studies with greater precision. However, many eC02 
studies do not report measures of variation in the data (standard error, standard deviation, or variance). 
In lieu of the measures of variance, studies can be weighted by the number of replicates (n) or, alter- 
natively, each study can be assigned equal weight, that is, unweighted method {Jablonski et al., 
2002). I used both methods (weighted and unweighted) to calculate the means of effect sizes with 
95% CIs and compared the results of both methods. Nearly in all instances, the difference between the 
weighted and the unweighted means was small and lesser than the standard error of the unweighted 
mean. For example, across all the FACE studies, the overall mineral change was -6.1% (-7.8 to -4.4) 
when unweighted cf. the -6.5% (-8.0 to -5.1) when weighted. For the reason of close similarity 
between weighted and unweighted approaches, I used the simpler out of the two methods, that is the 
unweighted one, when reporting the results. 

Since the distribution of effect sizes is not necessarily normal, I applied both parametric (ttest) and 
non-parametric (bootstrapping with 10,000 replacements) tests for calculating the 95% CI for the 
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mean effect size and the statistical power. The latter was calculated for: (1) an absolute effect size of 
5%, and (2) the probability of Type I error, a = 0.05. If the variance of a small sample << the true pop- 
ulation variance, then this leads to substantial overestimations of Cohen's d and the statistical power. 
To be conservative when estimating power for small samples (m <20), I used the larger of the sample 
standard deviation or 0.21 , which is the standard deviation for the entire mineral dataset. 

The results from the parametric and non-parametric tests were very close. For example, for Zn in 
edible tissues (sample size = 65), ttest yields (-11.4%, -14.0 to -8.7) and 0.91 power cf. (-11.4%, -13.9 
to -8.7) and 0.92 power for the bootstrapping procedure. A close similarity between the results of ttest 
and non-parametric test is expected when sample size (m, the number of independent observations for 
each mineral) is >30, which often was the case in this study. For reporting purposes, I used the 95% CI 
and the power generated by the non-parametric method, that is, the bootstrapping procedure. 

Testing for publication bias 

To test for publication bias or 'the file drawer effect' in the dataset, I plotted effect sizes against corre- 
sponding sample sizes/replicates, n, to provide a simple visual evaluation of the distribution of effect 
sizes {Figure 9). The resulting cloud of points is funnel-shaped, narrowing toward larger sample sizes, 
and overall is symmetrical along the mean effect size. This indicates the absence of any significant 
publication bias {Egger et a/., 1997). 

Fragmenting the dataset into categories 

Meta-analytic CO2 studies often partition their datasets into various categories (e.g., plant group, 
plant tissue, fertilization, or water regime) to estimate effect sizes for each category. Such data 
fragmentation, however, is warranted only if the statistical power of the resulting test for each 
category is adequate. Otherwise, low power can lead to non-significant outcomes and Type II 
errors. As tempting as it can be to partition the current dataset into many categories and cases 
(e.g., Zn in fruits, Fe in tuber, Cu in annuals, multiple CO2 levels), only by fragmenting the data into 
sufficiently large categories an adequate statistical power can be retained. Such categories 
include: foliar tissues, edible tissues, woody plants (trees and shrubs), herbaceous plants, FACE 
studies, non-FACE studies, crops, wild plants (all non-crops, including ornamental plants), C3 
plants, C4 plants, rice, wheat, barley, and potato. Furthermore, I fragmented the data for C3 plants, 
the foliar and the edible tissues, the non-FACE and the FACE studies into individual chemical ele- 
ments and into individual common plant names (e.g., all rice cultivars grouped under 'rice'). For 
the regional analysis, I used only OTC and FACE studies because they reflect local environment 
much more accurately than studies using complete-enclosures (e.g., closed chamber, glasshouse). 
If an OTC or FACE study did not report precise geographic coordinates, then the latitude and longi- 
tude of a nearby research facility or city was used (all coordinates in the dataset are in decimal units). 
Figures 1-7 include results with the statistical power >0.40 for each element, country, region, plant 
tissue or category. Generally, power >0.80 is considered acceptable (Cohen, 1988). Unfortunately, 
such a level was achievable only for elements for which the data are most abundant and for the 
ionomes of some plant groups and species. Note that the power was calculated for a 5% effect size, 
while the true effect size is likely to be larger (~8%); therefore, the true power is likely to be higher than 
the calculated power for most results. All the results, irrespective of the statistical power, can be found 
in Figure 1-source data 1. Furthermore, Figure 1 shows the mean effect sizes (with their 95% CI) 
plotted against their respective statistical powers for all the minerals and all the plant groups/tissues. 

Formula for calculating percentage changes in TNC:protein and 
TNC:minerals 

If the concentration of substance X in a plant increases by x% and concomitantly the concentration of 
substance Y decreases by y% in the plant, then the X-to-Y ratio of the plant (X:Y) increases by: 

X + V 

^ -100% (1) 



100-y 

Proof 

Let us denote the initial concentrations of substances X and Yin a plant as and y^, respectively. Suppose 
the Xand Ycontents in the plant changed by x% and -y%, respectively. Then the new X content in the 
plant, Xg, is 
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xe = xa-(100 + x)%, 

and the new / content In the plant,yE, is 

ye = yA "(1 00 -y)%. 

The original X : Y = xa/ya, while the new X : Y = xe/ye- Since the percentage change in the X:Yequals to: 

new 



new - original ^^^^^^^ 
original 



original 



1 



100%, 



substituting xaIya and xeIye for the original and the new, respectively, yields: 

Xe/ye _ ^ ^ XE ■ yE _ ^ ^ (xa (100 + x) %)■ yA _ ^ ^ 100 + x _^ ^ x + y 
xa/ya xa-Ya XA(yA(100-y)%) 100-y 100-y' 

An advantage of Equation 1 is that it holds true irrespective of whether the decrease in Y is driven 
by some reason applicable only to Y or by the increase in X, that is the dilution by X. 
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